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Abstract 

Recent experiments suggest that certain data of B — ► tttt, ttK decays are inconsistent 
with the standard model expectations. We try to explain the discrepancies with R-parity 
violating supper symmetry. By employing the QCD factorization approach, we study these 
decays in the minimal supersymmetric standard model with R-parity violation. We show 
that R-parity violation can resolve the discrepancies in both B — > tttt and B — > irK decays, 
and find that in some regions of parameter spaces all these requirements, including the CP 
averaged branching ratios and the direct CP asymmetries, can be satisfied. Furthermore, 
we have derived stringent bounds on relevant R-parity violating couplings from the latest 
experimental data, and some of these constraints are stronger than the existing bounds. 
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1 Introduction 



The detailed study of B meson decays plays an essential role for understanding the CP violation 
and the physics of flavor. Recent experimental measurements [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12] 
have shown that some hadronic B decays to two pseudoscalar mesons deviate from the standard 
model (SM) expectations . 

In the B — ► txtx decays, there are three such discrepancies: the direct CP asymmetry for the 
mode B — > tt + it~ is very large [2, 10], the 7r°7r° mode is found to have a much larger branching 
ratio (» 1.5 x 10~ 6 ) [4, 11] than the SM expectations (~ 10~ 7 ) [13, 14], and the theoretical 
estimation of B — > 7r + 7r~ branching ratio [13, 14] are about 2 times larger than the current 
experimental average [1, 7]. The il B — > irir puzzle" is reflected by the following quantities [15]: 
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here we use tb u /tb c1 = 1.069 [16], the central values calculated within the QCD factorization 
(QCDF) give = 1.24 and R$ = 0.07 [13]. In Ref. [15], Buras et al. pointed out these data 
would indicate the large nonfactorizable contributions rather than new physics (NP) effects, 
and could be perfectly accommodated in the SM. However, it is hard to be realized by explicit 
theoretical calculations. 

The B — > txK system consists of the four decay modes B d — > tx^K^^B^ — > tt^K^B^ — > 
^K^ 1 and B d — > 7r°K°, which are governed by QCD penguin process b — > sqq. The BABAR 
and Belle collaborations have measured the following ratios of the CP averaged branching ratios 
[17]: 
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where we have included the latest Belle [6] and BABAR [12] measurements. The ratio R, which 
is expected to be only marginally affected by color-suppressed electroweak (EW) penguins, does 
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not show any anomalous behavior. The "B — > 7rif puzzle" is reflected by the small value of 
R n which is significantly lower than R c . Since R c and Rn could be affected significantly by 
color-allowed EW penguins, the U B — > 7rfT puzzle" may be a manifestation of NP in the EW 
penguin sectors [15, 18, 19], and will offer an attractive avenue for physics beyond the SM to 
enter the B — > nK system [20]. 

Although these measurements represent quite a challenge for theory, the SM is in no way 
ruled out yet since there are many theoretical uncertainties in low energy QCD. The recent 
theoretical results for B — > irK [21, 22] show that the next to leading order corrections may be 
important. However, it will be under considerable strain if the experimental data persist for 
a long time. Existence of NP as possible solutions have been discussed in Refs. [15, 18, 23]. 
Among those NP models that survived EW data, one of the most respectable options is It- 
parity violating (RPV) supersymmetry (SUSY). The possible appearance of RPV couplings 
[24], which will violate the lepton and baryon number conservation, has gained full attention 
in searching for SUSY [25, 26]. The effect of RPV SUSY on B decays have been extensively 
investigated previously in the literatures [27, 28, 29, 30, 31, 32]. In this work, we extend our 
previous study of the B — > VV decays [31] to the B — > nn, ttK decays using RPV SUSY 
theories by employing the QCD factorization (QCDF) approach[33] for hadronic dynamics. At 
the same theoretical ground, it would be interesting to know if we could find solutions to the 
B — » 7T7T and B — > irK puzzles besides the B — > VV polarization puzzle [31]. We show that the 
B — » 7T7T, nK puzzles could be resolved in the presence of the RPV couplings. Moreover, using 
the latest experimental data and theoretical parameters, we try to explain all available data 
including the CP averaged branching ratios and the direct CP asymmetries by the relevant 
RPV couplings. We note that the branching ratios of the B — > tttt, ttK decays have also been 
studied in RPV SUSY in [32, 34]. In this study, we present a new calculation of the decays 
with up-to-date inputs, such as form factors, experimental measurements and so on. 

The paper is arranged as follows. In Sec. II, we calculate the branching ratios and the CP 
asymmetries of B — > tttt, ttK decays, which contain the SM contributions and the RPV effects 
using the QCDF approach. In Sec. Ill, we tabulate theoretical inputs in our analysis. Section 
IV is to deal with data and discussions, we also display the allowed regions of the parameter 
space which satisfy all the experimental data. Section V contains our summary and conclusion. 
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2 The theoretical frame for B — > tttt and 7ri^ decays 



2.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the AB = 1 transition at the scale \i is 
given by [35] 



G { 10 

^■eff = ~?K E K] ClQl + C 2 Q 2 + E \piQi + + CsgQt 

V ^ p=«,c I j=3 



•+/I.C, (1) 



here A p = V p bV* (p G {it, c}) and the detailed definition of the operator base can be found in 
[35]. 

The decay amplitude for — > MxM 2 is 

A SM (B -> MiM 2 ) = (MiM 2 |Wf^|5 

= %EE A P C^)( M i M 2|g^)|5). (2) 

The essential theoretical difficulty for obtaining the decay amplitude arise from the evaluation of 
hadronic matrix elements (MiM 2 \Qi{[i)\B) . There are at least three approaches with different 
considerations to tackle the said difficulty: the naive factorization (NF) [36, 37], the perturbative 
QCD [14], and the QCDF [33]. 

The QCDF [33] developed by Beneke, Buchalla, Neubert and Sachrajda is a powerful frame- 
work for studying charmless B decays. In Refs. [13, 38, 39], B — > 7nr,nK decays have been 
analyzed in detail in the SM with the QCDF approach. We will also employe the QCDF 
approach in this paper. 

The QCDF [33] allows us to compute the nonfactorizable corrections to the hadronic matrix 
elements {MxM2\Oi\B) in the heavy quark limit. The decay amplitude has the form 

A SM (B - M,M 2 ) = ^EE A p{ a n M 2|J 2 |0)(M 1 |J 1 | J B) +b p i {M 1 M 2 \J 2 \0)(0\J 1 \B)j, (3) 

where the effective parameters af including nonfactorizable corrections at order of a s . They 
are calculated from the vertex corrections, the hard spectator scattering, and the QCD penguin 
contributions. The parameters &f are calculated from the weak annihilation contributions. 

Following Beneke and Neubert [13], coefficients af can be split into two parts: af = af j+af u . 
The first part contains the NF contribution and the sum of nonfactorizable vertex and penguin 




Figure 1: The next to leading order nonfactorizable contributions to the coefficients a 1 -. 



corrections, while the second one arises from the hard spectator scattering. The coefficients 
read [13] 
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9n N c 

where a s = a s (fi), Cf = (N c — l)/(2Nc), and Nc = 3 is the number of colors. The quantities 
Vm, Hm 2 m x > Pm2i Pm,& Pm^ W an d Pjaz^ cons ist of convolutions of hard-scattering kernels with 
meson distribution amplitudes. Specifically, the terms Vm come from the vertex corrections in 
Fig. l(a)-l(d), P M 2 and P M3 (Pm^ W an d Pjfi ) ar ise from QCD (EW) penguin contractions 
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and the contributions from the dipole operators as depicted by Fig. 1(e) and 1(f). Hm 2 m 1 is 
due to the hard spectator scattering as Fig. 1(g) and 1(h). For the penguin terms, the subscript 
2 and 3 indicate the twist 2 and 3 distribution amplitudes of light mesons, respectively. 
In Eq.(4), Vm (M = ir,K) contains the nonfactorizable vertex corrections , which is 
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Next, the penguin contributions at the twist- 2 are described by the functions 
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where rif = 5 is the number of quark flavors, and s u = 0, s c = {m c jm^) 2 are mass ratios involved 
in the evaluation of the penguin diagrams. The function Gm(s) is defined as 
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The twist-3 terms from the penguin diagrams are given by 
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with 
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Finally, the hard spectator interactions can be written as 
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Considering the off-shellness of the gluon in hard scattering kernel, it is natural to associate a 



scale Uh ~ \ ^QCD m b > rather than /i ~ mj,. For the logarithmically divergent integral, we will 



parameterize it as in [13]: Xh = J^du/u = —Iti(Aqc d / 'mb) + Qh^ h mb/^QCD with {qh,<I>h) 
related to the contributions from hard spectator scattering. In the later numerical analysis, we 
shall take Aq C d = 0.5GeV, 4>h) = (0, 0) as our default values. 




Figure 2: The weak annihilation contributions to the coefficients 6f . 



At leading order in a s , the annihilation contribution can be calculated from the diagrams 
in Fig. 2. The annihilation coefficients (61,62), (63,64) and (bf w ,bf w ) correspond to the contri- 
butions of the tree, QCD penguins and EW penguins operators insertions, respectively. Using 
the asymptotic light cone distribution amplitudes of the mesons, and assuming SU(3) flavor 
symmetry, they can be expressed as 
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and 
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Here the superscripts i and f refer to gluon emission from the initial and final state quarks, 
respectively. The subscript k of A l £ refers to one of the three possible Dirac structures Ti ® T 2 , 
namely k = 1 for (V -A) ® (V - A), k = 2 for (V — A)® (V + A), and k = 3 for (-2)(S-P)<8> 
(S + P). Xa = Jo du/u is a logarithmically divergent integral, and will be phenomenologically 
parameterized in the calculation as Xh- As for the hard spectator terms, we will evaluate the 
various quantities in Eq.(12) at the scale /x^ = J KqcD^b- 



With the coefficients in Eq.(4) and (11), we can obtain the decay amplitudes of the SM part 
Af AI (the subscript "/" denotes the part of QCDF contribution) and A^ M (the subscript "a" 
denotes the annihilation part). B — > titt and irK decay amplitudes are given in Appendix A. 

2.2 R-parity violating SUSY effects in the decays 

In the most general superpotential of the minimal supersymmetric Standard Model (MSSM), 
the RPV superpotential is given by [40] 

V\fc = mLiH u + + KikUQjDt + l% w UiDjBk, (13) 

where L and Q are the SU(2)-doublet lepton and quark superfields and E c , U c and D c are the 
singlet superfields, while i, j and k are generation indices and c denotes a charge conjugate field. 

The bilinear RPV superpotential terms [iiLiH u can be rotated away by suitable redefining 
the lepton and Higgs superfields [41]. However, the rotation will generate a soft SUSY breaking 
bilinear term which would affect our calculation through penguin level. However, the processes 
discussed in this paper could be induced by tree-level RPV couplings, so that we would neglect 
sub-leading RPV penguin contributions in this study. 

The A and A' couplings in Eq.(13) break the lepton number, while the A" couplings break 
the baryon number. There are 27 A^- fc couplings, 9 A^ and 9 A-' jfc couplings. X\%j]k are antisym- 
metric with respect to their first two indices, and A^ fe i are antisymmetric with j and k. The 
antisymmetry of the baryon number violating couplings A|L fe i in the last two indices implies 
that there are no X"j k operator generating the b — > sss and b — * ddd transitions. 

u i u m dj d m 

X. e-u / \. v u 

/ A ijk \ mn \ / A ijk 

dk d n ^ k d n 

Figure 3: Sleptons exchanging diagrams for nonleptonic B decays. 

From Eq.(13), we can obtain the following four fermion effective Hamiltonian due to the 
sleptons exchange as shown in Fig. 3 
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Figure 4: Squarks exchanging diagrams for nonleptonic B decays. 
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The four fermion effective Hamiltonian due to the squarks exchanging as shown in Fig. 4 are 
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and Pq = 11 — gTi/. The subscript for the currents 
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(jn)i,8 represents the current in the color singlet and octet, respectively. The coefficients rj 
and r]~ 8 ^ are due to the running from the sfermion mass scale rrij (100 GeV assumed) down 
to the m;, scale. Since it is always assumed in phenomenology for numerical display that only 
one sfermion contributes at one time, we neglect the mixing between the operators when we 
use the renormalization group equation (RGE) to run down to the low scale. 

Generally, the RPV couplings can be complex and their phases may induce new contribu- 
tions to the CP violation, so we write them as 
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here is the RPV weak phase, which could be any value between and ir. To include the effect 
of 7r < 0£ < 27r, |Ay/-A* mn | is allowed to take both positive and negative values for simplicity. 

the Tfl . For B — > PP decays, since 
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(P| q lfl (l - l5 )b I B) = (P| g 7|t (l + l5 )b I B) = (P| q^b | B), 



(17) 
(18) 
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the RPV contribution to the decay amplitude will modify the SM amplitude by an overall 
relation. 

Since we are considering the leading effects of RPV, we need only evaluate the nonfac- 
torizable vertex corrections and hard spectator scattering contributions. We ignore the RPV 
penguin contributions, which are expected to be small even compared to the SM penguin am- 
plitudes, due to the smallness of the relevant RPV couplings compared with the SM gauge 
couplings. As shown in Ref. [34], the bounds on the RPV couplings are insensitive to the 
inclusion of RPV penguins. We also have neglected the annihilation contributions in the RPV 
amplitudes. The R-parity violating part of the decay amplitudes can be found in Appendix 
B. 

2.3 The branching ratio and the direct CP asymmetries 

With the QCDF approach, we can get the total decay amplitude 

A{B -> P 1 P 2 ) = A s f M {B -> P 1 P 2 ) + A s a M {B -> P 1 P 2 ) + A${B^ P 1 P 2 ). (19) 

The expressions for the SM amplitude A s ^ and the RPV amplitude A? are presented in 
Appendices A and B, respectively. From the amplitude in Eq. (19), the branching ratio reads 

Br(B - P 1 P 2 ) = ^4 \A(B - P.P^fS, (20) 



iirm 



where S = 1/2 for identical P\ and P 2 , S = 1 otherwise, tb is the B lifetime, \p c \ is the center 
of mass momentum of light mesons in the rest frame of B meson, and given by 

\Pc\ = 2~\/l m B ~ ( mp i + m Pi) 2 ][ m B - ( m Pi - m P 2 ) 2 ]- ( 21 ) 

The CP averaged branching ratios are defined by 

Br(B ± -> P±P 2 ) = - [Br(B~ -> PfP 2 °) + Br{B + -> P?P%)] , 
Br{B d ->P±P?) = l -\Br{B' d ^P+P 2 )+Br{B Q d ^P{P+)\, (22) 



Br(B d ^P?P°) = -lBr(B d ^P?P°)+Br(B d ^P?P°) 

The direct CP asymmetry is defined by 

dir _ Br(B -> f) - Br(B -> f) 

~ Br(B ->•/)+ Br(B /) ' {Z6) 
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3 Input parameters 



A. Wilson coefficients 

For numerical analyzes, we use the next-to-leading Wilson coefficients calculated in the naive 
dimensional regularization (NDR) scheme and at mj scale 



C x = 1.082, C 2 = -0.185, C 3 = 0.014, C 4 = -0.035, C 5 = 0.009, 
C 6 = -0.041, C 7 /a e = -0.002, C 8 /a e = 0.054, C 9 /a e = -1.292, 
Cio/ae = 0.263, C$ = -0.299, C& f = -0.143. 



B. The CKM matrix element 

The magnitude of the CKM elements are taken from [42] 



\Vub\ 

\v cb \ 



0.00367 ±0.00047, 



0.0413 ±0.0015, 



(24) 



\V ud \ = 0.9738 ± 0.0005, \V US \ = 0.2200 ± 0.0026, 
\Ved\ = -0.224 ± 0.012, IV^I = 0.996 ± 0.013, 
\V t * b V td \ = 0.0083 ± 0.0016 \V tb V t * s \ = -0.047 ± 0.008, 

and the weak phases 7 = 60° ± 14°, sin(2/3) = 0.736 ± 0.049. 
C. Masses and lifetime 

There are two types of quark mass in our analysis. One type is the pole mass which appears 
in the loop integration. Here we fix them as 



m u = m d = m s 



0. 



1.47 GeV, m b = 4.8 GeV. 



The other type quark mass appears in the hadronic matrix elements and the chirally enhanced 
factor = =^ through the equations of motion, which are renormalization scale dependent. 
We shall use the 2004 Particle Data Group values [42] for discussion (the central values are 
taken as our default values) 

m u (2GeV) = 0.0015 ~ 0.004 GeV, m d (2GeV) = 0.004 ~ 0.008 GeV, 
m s (2GeV) = 0.08 ~ 0.13 GeV, m b (m b ) = 4.1 ~ 4.4 GeV, 



and then employ the formula in Ref. [35] 
m(n) = m(/i ) 
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to obtain the current quark masses at /i scale. The definitions of jffl ,y$ , Po, Pi can be found 
in [35]. 

D. The light cone distribution amplitudes (LCDAs) of the pseudoscalar meson 

For the LCDAs of the pseudoscalar meson, we use the asymptotic form [43, 44] 

$ P (x) = 6x(l - x), $p fx) = 1. (25) 

We adopt the moments of the $f (£) defined in Ref. [13, 33] for our numerical evaluation 

'd^=^, (26) 

with Xb = 0.46 GeV [45]. The quantity A# parameterizes our ignorance about the B meson 
distribution amplitudes and thus brings considerable theoretical uncertainty. 

E. The decay constants and form factors 

For the decay constants, we take the latest light-cone QCD sum rule (LCSR) results [46] in our 
calculations 

f B = 0.162 GeV, f K — 0.16 GeV, f K = 0.131 GeV. 

For the form factors involving B — > K and B — > 7r transitions, we adopt the center values of 
the results [46] 

f^ K (0) = 0.331, fo^ n (0) = 0.258. 



4 Numerical results and analysis 

We will present our numerical results in this section. At first, we will show our estimations 
in the SM by taking the center value of the input parameters and compare with the relevant 
experimental data. Then, we will consider the RPV effects and constrain the relevant RPV 
couplings by the averages of Belle and BABAR measurements of the CP averaged branching 
ratios and the direct CP asymmetries. 

When considering the RPV effects, we will use the input parameters and the experimental 
data which are varied randomly within la level and 2a level, respectively. In the SM, the weak 
phase 7 is well constrained, however, with the presence of RPV, this constraint may be relaxed. 
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We would not take 7 within the SM range, but vary it randomly in the range of to tt to obtain 
conservative limits on RPV couplings. We assume that only one sfermion contributes at one 
time with a mass of 100 GeV. So for other values of the sfermion masses, the bounds on the 
couplings in this paper can be easily obtained by scaling them by factor f 2 = ( \0QGeV ) 2- 

The main numerical results in the SM and the relevant data from the Belle collaborations 
[7, 8, 9, 10, 11, 12] and BABAR collaborations [1, 2, 3, 4, 5, 6] are presented in Table I, II and 
III, which show the results for the CP averaged branching ratios, the direct CP asymmetries 
and the ratios of the CP averaged branching ratios, respectively. 

From Table I, II, III, we can see the puzzle in B — > tttt, ttK which we have already men- 
tioned in the introduction. For example, the new experimental data for B — > 7r°7r°, 7r 
branching ratios are significantly larger than the SM predictions, moreover, the expected rela- 
tion Acp(B° — > 7T ± A' :F ) ps A^jpiB^ — > 7r°A' ± ) obviously contradict to the experimental data, 
even with the opposite sign for them, the value of R n is significantly lower than R c , and so on. 

Now we turn to the RPV effects which may give possible solutions to the puzzle. We use 
the CP averaged branching ratios, the direct CP asymmetries and the relevant experimental 

Table I: The CP averaged branching ratios of B — > ix-k^hK in the SM(in unit of 10 -6 ). 
Experimental data from BABAR and Belle and the SM predictions in the framework of NF 
and QCDF, where Br^ and Br^ +a denote the results without and with the contributions 
from weak annihilation, respectively. 



Decays 


Belle 


BABAR 


Average 


NF 


QCDF 


Br f 




B d - 


■> 7T ± 7T =F 


4.4 ±0.6 ±0.3 


4.7 ±0.6 ±0.2 


4.56 ±0.46 


7.56 


7.88 


8.33 


B ± - 


^ttV 


5.0 ±1.2 ±0.5 


5.8 ±0.6 ±0.4 


5.61 ±0.63 


5.19 


4.90 


4.90 


B d - 


^7T 7T 


9 0+0.4+0.2 

z -°-0.5-0.3 


1.17±0.32±0.10 


1.45 ±0.29 


0.17 


0.16 


0.18 


B ± - 


-> 7T ± K 


22.0 ± 1.9 ± 1.1 


26.0 ± 1.3 ± 1.0 


24.57 ± 1.31 


12.42 


14.77 


17.09 


B ± - 




12.0 ± 1.3±J;jj 


12.0 ±0.7 ±0.6 


12.00 ±0.81 


7.23 


8.36 


9.57 


B d - 


* 7T ± A' T 


18.5 ± 1.0 ±0.9 


17.9 ±0.9 ±0.7 


18.15 ±0.87 


9.93 


11.64 


13.54 


B d - 


+ 7l°K 


11.7d= 2.3±i;| 


11.4 ±0.9 ±0.6 


11.44 ± 1.00 


4.20 


5.06 


5.97 
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Table II: The direct CP asymmetries A^p (in unit of 10~ 2 ) for B — > irn, ttK. 
Experimental data from BABAR and Belle. A.q P and are defined in the 

similar way as branching ratios in Table I. 



Decays Belle BABAR Average 



QCDF 



\f Af+a 
•^CP -^CP 



Bd^Tr^ 58 ±15 ±7 9 ±15 ±4 31.93 ± 11.32 -5.86 -5.68 

B ±^ n o n ± 2 ± 8 ± 1 -1±10±2 0.85 ±6.32 -0.07 -0.07 

B d -> 7T 7r° 44i| ±17 12 ± 56 ± 6 28.33 ± 39.42 63.04 60.55 

B ± ^n ± K 5±5±1 -9 ± 5 ± 1 -2.00 ± 3.61 1.23 1.13 

B ± ^7T°K ± 4±4±2 6 ± 6 ± 1 4.70 ± 3.60 8.12 7.29 

B d ^n ± KT -11.3 ±2.2 ±8 -13.3 ±3± 0.9 -12.02 ± 1.86 6.32 5.50 

B d ^n K -12 ±20 ±9 -6 ±18 ±6 -8.57 ± 14.35 -2.45 -2.19 



Table III: The ratios of the CP averaged branching ratios. 



Ratios 


Exp. 


NF 




QCDF 


without ann. with ann. 




2.30 ±0.35 


1.28 


1.16 


1.10 


D7T7T 

-"-oo 


0.66 ±0.14 


0.04 


0.04 


0.04 


R 


0.79 ±0.06 


0.80 


0.79 


0.79 


R c 


0.98 ±0.08 


1.16 


1.13 


1.12 


Rn 


0.79 ±0.08 


1.18 


1.15 


1.13 



averages of Belle and BABAR to constrain the spaces of the RPV parameters. As known, data 
on low energy processes can be used to impose rather strictly constraints on many of these 
couplings. The random variation of the parameters subjecting to the constraints as discussed 
above leads to the scatter plots displayed in Fig. 5 and 6. 

The B — > 7T7T decays involve the quark level processes b — > dqq (q = u,d). The all three RPV 
couplings maybe resolve the B — > titt puzzle, which has been shown in Fig.5. These allowed 
parameter spaces are essentially controlled by the CP averaged branching ratios and the direct 
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-1.0 -0.5 0.0 0.5 1.0 



h, 3 O x1 ° 3 

Figure 5: The allowed parameter spaces for the relevant RPV couplings constrained by B — ► tttt. 

CP asymmetries of different tttt modes. From Fig. 5, we can see that RPV weak phase is not 
constrained so much, but the relevant magnitude of the RPV couplings are constrained within 
rather narrow ranges. We can obtain the allowed parameter spaces for the relevant couplings, 
which are summarized in Table IV. For comparison, we also list the existing bounds on these 
quadric coupling products [25, 32]. 



Table IV: Bounds for the relevant RPV coupling constants by B — > tttt decays 
for 100 GeV sfermions and previous bounds are listed for comparison. 



Couplings 


Bounds 


Process 


Previous bounds 


l A 132 A 112l 


[5.4 x 


10" 4 , 2.9 x lO" 3 ] 


B d - 


-> 7r ± 7r =F , 7r°7r° 


< 5.0 x 10~ 3 [25] 


I\l3-Mlll 


[2.7 x 


10" 4 , 7.7 x 10" 4 ] 


B d - 


-> 7l^~7T~^ , 7r 7T 


< 2.5 x 10~ 3 [25] 


1-^11^31 1 


[3.8 x 


10" 4 , 1.0 x 10~ 3 ] 


B - 




< 2.0 x 10" 3 [32] 



We note that since the quark content of n° is antisymmetric combination (uu — dd) / a/2, 
the decays B — > n^ii^, 7r°7r° could be induced by superpartners of both up-type and down- 
type fermions in RPV SUSY. For example, b — > ddd could be induced by sneutrino, while 
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— — \f \/* 

b — > uud could be induced by slepton with the same 

X' ilz \'*u product. We take ' 13 ? n and 
il3 . 2 ni contribute to b — >• uud and b — > ddd at the same time, and the effects of A'ioA'11 will 

771- " 111 

be summed. So that A-^A-^ in amplitudes of B ± — > 7T°7r ± is cancelled and in amplitudes of 
B c i — > 7t°it is partly cancelled if taking = mo i . 

The B — > 7rA' processes are due to 6 — > sqq (q = w, d) transitions at quark level. There 
are six RPV coupling constants contributing the four B — > txK decay modes. We scan the 
parameter space for possible solutions, and find that only four pairs of RPV couplings, A" 31 Af 21 , 
^131^121' -Mi3-Mi2 an d A^iA^ 31 , can survive after satisfying all relevant experimental data of 
B — > nK decays. But we do not get the solutions to the experimental data at 2a level with the 
other two pairs of RPV coupling constants, A^ 23 A^ n and A£ n A£ 32 . Figure 6 displays the allowed 
ranges for RPV couplings which satisfy all relevant experimental data of the B — > irK decays. 
The constraints for the four RPV couplings are summarized in Table V. For comparison, we 
also list the pervious bounds [25, 32]. 




Figure 6: The allowed parameter spaces for the relevant RPV couplings bounded by B — > ttK. We 
have not listed the RPV couplings, A^ 23 A- n and A£ n A^ 32 , for which we have not found the allowed 
space. 
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Table V: Bounds for the relevant RPV coupling couplings for 100 GeV sfermions by 
B — » irK decays and previous bounds are also listed. 



Couplings 


Bounds 


Process 


Previous bounds 


1 X" X" 1 

l A i31 A i2ll 


[13.6xl(T 3 , 15.8X10" 3 ] 
[3.7xlCT 4 , 4.6xlCT 3 ] 


B ± -> tt ± K, 7r° J fs: ± 


< 4. x 10" 3 [25] 


1 X" X" 1 
1 A 131 A 12ll 


[5.6xlCT 3 , 7.2xl(T 3 ] 
[6.IXIO- 4 , 4.6xl0~ 3 ] 


B -> K°K ± ,K ± K*,n°K 


< 4. x 10~ 3 [25] 


IX' X' 1 

1 A il3 A il2l 


[3.6 x 10~ 4 , 1.1 x 10~ 3 ] 


four modes of 7rK 


< 5.7 x 10~ 3 [25] 


IX' X' 1 

l A i21 A i3ll 


[l.lxlO -3 , 1.2xl0~ 3 ] 
[6.IXIO- 4 , 8.7xl0" 4 ] 


B -> 7T ± K, TT ^, 7T°K 


< 1.3 x 10~ 3 [32] 



For the coupling A" 31 A" 21 , we get two ranges by the CP averaged branching ratios and the 
direct CP asymmetries of B — > irK. The bound for this coupling is 10" 1 ~ 10~ 2 by branching 
ratios of B -> irK in [32], and the limit of Af 3fc A'/ 2fc is 0.16 by B -> K*^ decay [25]. So both 
spaces of A" 31 A" 21 may be allowed. The A'/^A'/^ also have two ranges, we may get only the space 
1^131-^121 1 e [0-61 x 10~ 3 ,4.60 x 10~ 3 ] if considering the constraint for A" 31 A" 21 . The bounds for 
couplings A^ 13 A^ 12 and A^ 21 A^ 31 are obtained by four or three decay modes of B — > irK, their 
ranges are also very narrow. For the RPV couplings A- 23 A- n and A- n A- 32 , we have not obtained 
their solutions to the puzzles. 

The B — > ir°K ± and B° — > tt°K° decays could be induced by superpartners of both up- 
type and down-type fermions, b — > dds could be induced by sneutrino, while b — > uus could be 
induced by slepton with the same A^ 13 A-^ 2 product. For the same reason as in B — > tt ^, 7r 7r°, 
the effects of A^ 13 A^ 2 have been summed. 

The above analysis has shown that the puzzles in the B — > tttt, irK decays can be resolved 
with RPV effects, however, the solution parameter spaces are always very narrow. The allowed 
spaces constrained by the B — > tttt, ttK decays are consistent with that by B — > VV decays in 
our previous study [31]. 

5 Conclusions 

The recent observations of B — > tttt, ttK decays which are inconsistent with the SM expectations 
represent a challenge for theoretical interpreting. We have employed the QCDF to present a 
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study of the RPV effects in the B — > 7T7t, txK decays. Our analysis has shown that a set 
of RPV couplings play an important role to resolve the discrepancies between the theoretical 
predictions in the SM and the experimental data. However, the windows of the RPV couplings 
intervals are found to be always very narrow. It implies that these couplings, part or all of 
them, might be pinned down from the rich experimental phenomena in these decays. However, 
it also implies the window could be closed easily with refined measurements from experiments 
in the near future. 

It should be noted that some of the A' couplings can generate sizable neutrino masses [41, 47]. 
Allanach et al. have obtained quite strong upper bound A^ < 10~ 3 in the RPV mSUGRA 
model. Our bounds for quadric products |A^ 13 A- n | and |A- 11 A^ 31 | are of order of 1CT 3 — 10~ 4 . 
So combining their constraints from neutrino masses and ours from B — > hit decays, the A' 
resolution window might be closed. However, we note that the constraints on A' from neutrino 
masses would depend on the explicit neutrino masses models with trilinear couplings only, 
bilinear couplings only, or both [41]. 

Furthermore, the A' m coupling has been constrained as low as |A' m | < 5.2 x 10~ 4 by neutrino- 
less double beta decay [48]. There are also strong bounds \X" 21 \ < 10" 4 and \X" 31 \ < 10~ 4 from 
double nucleon decay [49] and neutron oscillations [49, 50], respectively. Combining these strong 
bounds, our solutions with A' 113 A' m , A' m A' 131 and A'/^A'^ RPV products should be excluded. 
However, from the comprehensive collation of bounds upon trilinear RPV couplings in Ref.[51], 
our other solutions still remain. Explicitly, the B — > -kit puzzle could be resolved by the 
presences of RPV couplings Xi 32 X'{ 12 , A^ 11 A^ 31 and A- 13 A- 11 with i ^ 1, while the B — > ixK puzzle 
could be resolved by A-' 31 A-' 21 (z ^ 1), A' a3 A' a2 and A- 21 A- 31 . 

Generally, we can believe that QCDF calculations for the direct CP asymmetries could 
be much more accurate than that for the branching ratios, since many uncertainties could be 
cancelled in the ratios. Therefore the constraints from the direct CP asymmetries would be 
more well-founded than those only from branching ratio measurements [32]. Comparing our 
prediction with the recent experimental data within 2a level about the CP averaged branching 
ratios and the direct CP asymmetries, we have obtained bounds on the relevant products of 
RPV couplings. With more data from BABAR and Belle, one can significantly shrink the 
allowed parameter spaces for RPV couplings. We find that these constraints are consistent 
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with the previous bounds, even most of them are stronger than the existing limits [25, 31, 32], 
which may be useful for further study of RPV phenomenology. 

To summarize, we have shown that the B — > txtx puzzle and the B — > txK puzzle could be 
resolved in the RPV SUSY. Using the latest experimental data, we get the allowed values of the 
relevant RPV couplings, and the most of these new constraints are stronger than the existing 
bounds. 
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Appendix 



A . The amplitudes in the SM 



The factorized matrix elements are defined by 



A MlM2 = (M 2 |(g 1 g 2 ) (y _A)|0)(M 1 |(g 3 &)^ (27) 
B MlM2 ee (M 1 M 2 \(q 1 q 2 ) (v _ A) \0}(0\(q 3 b) {v . A) \B}=tf B f Ml f M2 . (28) 



A 6 ; M (B° -> 7T+7T 



Gi 



—7= {A Mrf ai + \ pd a\ + a? + r£(a£ + ag) j A 



TX TX 



0\ 



x/2 
G p 



■ ^ud( a l + a 2) + 2^ pd 



-a 7 + a 9 + a^ + r^af 



(29) 
(30) 



A T M (B° -> TX°TX ) 



_Gp( 
s/2 I 



pd 



3 



A~ 
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Af\B- 
A s f M {B- 



A s f M (B° 



ir-K°) 



A s f M (B° — > tt + K~ 



71 



°K°) 



A SM {B _ 


-> 7T + 7T ) 


A s a M (B- - 


-> 7T-7T ) 


A SM {B 


- 7r°7r°) 


A™(B- 


> 71- K°) 


A s a M (B- - 




A s a M (B° - 




a sm {b o _ 


7T°i? ) 



f / \ r ( p * p \ 

—j= <yA ps ya 4 -a 10 j 
G 

~Y {A us ai + A ps (af + a?„) + (ag + ajj) } 
+-^- |A us a 2 + ^A ps (-a 7 + a 9 )| A^, 

G 

— | {A us ai + X ps {a\ + a? ) + ?"^(ag + a s) } Artf, 



2 |Aps 
G 



04 - -c% ) + r x I a£ - -a\ 



1 

— ( 
2 



.4 



^ A n d&l + A 



63 + 26 4 - -b 3 



V2\ 
0, 

-^J= ^udW + A p d 

^{X U sb2 + X P s(h + bi w )}B 7Tl< 

^-{\usb2 + \ P s(h + bf W )}B„ K 



2~ 3 ' 2' 



2 4 



B 



ttK 



B T: 



B-„ 



ttK, 



(31) 
(32) 

(33) 
(34) 

(35) 

(36) 
(37) 
(38) 

(39) 

(40) 
(41) 

(42) 



here we note that A 



ps 



V pb V* s and X pd = V pb V* d . 




8\/2m~ 16V2m~. 



d 

r]- 8 ^ L K7 ,A K% 



V* V V* V 

A il3 A il2 A il3 A il2 



' \l* \l \l* \l 

A i2 , A A ill A ill A i32 



-8/Po n a 
1 1 i 



\8V2ml Lt 8v/2m? Li/ 



8\/2m~ T 8v2m~ r 

\//* \// W* W 

8mj 8m~ A 

\//* \// w* W 

A$(R° —i tt ^ - A 131 A 121 -4/A) p 4 A il3 A il2 -8/fl) r 4 

yl — »• 7r A J — — ——7= — ^-rj f K7T A Kn — — -= — — rj L K7T A KlT 

8v2mj 8v2m~ r 

d e Lt 

1* \r \t* \! \ 



v 8V2m? K 8V2m? M/ 
In the .4^, F MlM2 and A Mi m 2 are defined as 

F MlM2 = \- — + ^-9L{V M2 + H M2Mi \ 



Lm!M 2 = TT 



l ~^% {l2 + VKh + HhhM ^ 
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